Quantitative imaging of nanometric optical path length modulations by time-averaged 
heterodyne holography in coherent frequency- division multiplexing regime. 
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We report a demonstration of amplitude and phase imaging of out-of-plane sinusoidal vibration 
at nanometer scales with a heterodyne holographic interferometer. Time-averaged holograms of a 
phase-modulated optical field are recorded with an exposure time much longer than the modulation 
period. Optical heterodyning, a frequency- conversion process aimed at shifting a given radiofre- 
quency optical side band in the sensor bandwidth, is performed with an off-axis and frequency-shifted 
optical local oscillator. The originality of the proposed method is to make use of a multiplexed local 
oscillator to address several optical side bands into the temporal bandwidth of the sensor array. This 
process is called coherent frequency-division multiplexing. It enables simultaneous recording and 
pixel-to-pixel division of two side band holograms, which permits quantitative mapping of the mod- 
ulation depth of local optical path lengths yielding small optical phase modulations. Additionally, a 
linear frequency chirp ensures the retrieval of the local mechanical phase shift of the vibration with 
respect to the excitation signal, when screening a given frequency range. The proposed approach is 
validated by quantitative motion characterization of the lamellophone of a musical box, behaving as 
a group of harmonic oscillators, under weak sinusoidal excitation. Images of the vibration amplitude 
versus excitation frequency show the resonance of the nanometric flexural response of one individual 
cantilever, at which a phase hop is measured. 



Imaging nanometric optical path length modulations 
is of interest in non-destructive testing of micro electro- 
mechanical systems Among common approaches 
used for single-point nanometric vibration spectra 
measurements, laser Doppler heterodyne schemes Q 
are methods of choice, because of their high sensitivity 
and wide detection bandwidth, though imaging is 
hindered because scanning of the probe is required. 
In wide-field imaging of sinusoidal optical path length 
modulations, either performed with incoherent or 
coherent light sources, in direct image or holographic 
recording conditions, the mechanical phase is usually 
retrieved from phase-stepped stroboscopic schemes [2, [J. 
In the particular case of holographic measurements, 
this approach is extensively used in the sinusoidal 
excitation regime 0-[s|- Yet, in holography, the optical 
phase retardation of the recorded scattered field can be 
exploited for the determination of the local mechanical 
phase of an object motion with respect to the excitation 
signal when the readout rate of the array sensor is faster 
than the modulation period 0, [loj. An intermediate 
recording regime is quasi- time- averaging [3j, for which 
the modulation period is of the same order of magnitude 
as the frame rate. Nevertheless, because the vibration 
frequencies of interest in non-destructive vibration 
screening are still much higher than camera frame rates, 
time-averaged interferograms are recorded, even 
with high-speed cameras. Quantitative imaging of 
nanometer-scale optical path length modulations and 
phase retardation in stationary regimes, at the time 
scale of the frame exposure, may thus find applications. 

In this letter, we report a heterodyne holographic 
arrangement designed for amplitude and phase imag- 
ing of nanoscale out-of-plane sinusoidal vibrations in 
time- aver aged recording conditions without strobe light. 
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FIG. 1: Experimental set-up. 



The proposed arrangement consists of a Mach-Zehnder 
holographic interferometer in an off- axis and frequency- 
shifting configuration, with a multiplexed local oscillator 
(fig. [I]). Optical local oscillator multiplexing is realized 
by the coherent addition of radiofrequency (RF) signals 
driving an acousto-optic modulator in the reference 
beam channel. Coherent frequency-division multiplexing 
is a technique by which the total available bandwidth 
is divided into non-overlapping frequency sub-bands 
carrying separate signals. Addressing simultaneously 
distinct domains of the spectral bandwidth of holograms 
was performed with success for angular 11 5L wavelength 
[l6| . space-division [13], and side band [l8] multiplexing 
schemes. Coherent frequency-division multiplexing in 
heterodyne holography offers the opportunity to measure 
complex- valued maps of the RF spectral components of 
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the optical radiation field with accuracy and sensitivity, 
as long as they can be isolated in sub bands of the 
sensor's temporal bandwidth. In the reported exper- 
iment, high-quality factor flexural responses at audio 
frequencies of a lamellophone (the metallic comb of 
cantilevers from a musical box) excited in a sympathetic 
resonance regime are sought. Sympathetic vibration 
is a harmonic phenomenon wherein a vibratory body 
responds to external excitation in the vicinity of a 
resonance frequency. 

The optical setup used to perform non-destructive test- 
ing of the cantilevers' out-of-plane flexural responses at 
nanometric scales is sketched in figd It is designed to 
screen quantitatively steady-state vibration amplitudes 
and to detect mechanical phase jumps with respect to the 
excitation signal, by interferometric detection of an opti- 
cal field scattered by the lamellophone, referred to as the 
object field, beating against a reference optical field. The 
out-of-plane motion of the cantilevers provokes an optical 
path length modulation of the backscattered object field. 
Hence the temporal part of the object field E undergoes a 
sinusoidal phase modulation <j)(t) = 4nz(t)/\. The equa- 
tion of the out-of-plane motion z{t) of a cantilever, mod- 
eled as a damped harmonic oscillator of spring constant 
fc, mass m, damping coefficient c, driven sinusoidally by 
the force F(t) — Fo sm(ujt) provided by a piezo-electric 
actuator is 
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c^: + kz = F(t) 



(1) 



We seek a steady-state solution at the excitation fre- 
quency with an induced phase change of ip with respect 
to the excitation F(t), as a result of viscous damping 

z(t) = zo sin (cut + ip) (2) 

Its amplitude, zo, is proportional to the driving force 
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Two intermediate variables were introduced : the angular 
resonance frequency ujq of the undamped oscillator, and 
the dimensionless damping ratio £. They satisfy the re- 
lationships k = mujQ and c = 2muJo(. The phase shift of 
the sinusoidal response with respect to the driving force 
is 



ip = arctan 



(4) 



The resonance frequencies uj r = ujo yjl - 2C 2 of each can- 
tilever are almost equal to the resonant frequencies ujo 
of the undamped system. The optical phase modulation 
4>(t) = 0o sin(c<jt + ip) of depth </>o = 47rzo/A at the angu- 
lar frequency cj, can be decomposed on a basis of Bessel 
functions of the first kind J n (</>o)? via the Jacobi- Anger 
identity 



E = Se 
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FIG. 2: Magnitude of a statically-scattered light hologram 
\Ho\ (a). Magnitude of side band holograms (b) and 

| .Hi | (c). Amplitude map of the flexural response zo of the 
first cantilever for the first resonance at oj/(2tt) — 541 Hz 
(d). Phase images ip calculated from 128 raw interferograms 
around 540.9 Hz (e) and 541.8 Hz (f), in the neighborhood of 
the first resonance. Movies of the amplitude and phase maps 
of the out-of-plane vibration are reported in media 1 and 2. 



The quantity 



(6) 



is the complex weight of the optical side band of order 
n, where £ is the complex amplitude of the optical field, 
and (j)o = 47rzo / A is the modulation depth of the optical 
phase of the radiation backscattered by the cantilever. 
The local amplitude zo of the out-of-plane vibration and 
phase retardation ip with respect to the sinusoidal driv- 
ing force can be derived from the complex weight of two 
optical side bands with eq. [6] and the first-order Tay- 
lor developments of Jo(^o) and Ji(0o), near 0o = 0, for 
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The holographic interferometer is used for the detec- 
tion of RF phase modulation side bands of an object 
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FIG. 3: Magnitude of the FFT-spectrum (eqECQand e q. [T2l) . 
where three modulation side bands are addressed. 



field E in reflective geometry, beating against a LO field 
£?lo, as reported in [18|. Its purpose is to enable a ro- 
bust quantitative measurement of the nanometric mod- 
ulation amplitudes of cantilever resonances, as well as 
the associated phase shift with respect to the sinusoidal 
excitation signal. To do so, the adopted strategy is to 
separate the sensor bandwidth [— c^s/2, ojs/2[ into three 
non-overlapping frequency sub-bands, centered at 5uoi, 
t)co>2, and 6003, each of which is used to carry a separate 
optical modulation side band. The LO signal is composed 
of the sum of three phase-locked RF signals, to yield an 
optical LO field of the form 



E L o = W iu;Lt J>„ ( 



jAuJnt 



(8) 



where Au n = (n — 2)u — Su n is the frequency shift and 
c^n^LO is the complex weight of the LO component of 
rank n = 1,2,3. The positive parameters a n are the 
normalized weights of each LO component; in this ex- 
periment ai = 0L2 = as = 1/3. The LO components are 
tuned around the optical modulation bands of order —1, 
0, and +1. Precisely, the LO component of rank n beats 
against the optical modulation band of rank n — 2 at the 
apparent frequency Suj n , which is set within the camera 
bandwidth. Image rendering from the raw interferograms 
X(t) = \E(t) + £?lo(£)| 2 impinging on the sensor array 
involves a diffraction calculation performed with a nu- 
merical Fresnel transform. Its implementation in off-axis 
holographic conditions [HI yields complex- valued time- 
averaged 11 13 holograms I(t) from which the off- axis 
region 



H(t) = E(t)El Q (t) 



(9) 



is the support of the object-reference fields cross-term. 
The demodulation process (eq. [TT]) from N of consecutive 
recordings of finite exposure times acts as a band-pass 
filter of width us/N in the temporal frequency domain. 
Only the contributions within the sensor temporal band- 
width, between the Nyquist frequencies ±cjs/2, have to 
be taken into account in H(t), which can hence be ex- 
pressed as 



H(t) = ^2 £n-2 £l Q a n exp(iSuj n t) 



(10) 



where the quantities Suj n satisfy —uo$/2 < 5uo n < uo$/2. 
The FFT-spectrum of N consecutive off-axis holograms 

H K) = Ep=iH(27rp/oj s ) exp (-2inpir/N) (11) 

exhibits three peaks at the frequencies Slji, 8002, and Sus 
(fig. [3]). Their complex values yield holograms H m , of 
rank m = —1,0, 1, whose magnitude is proportional the 
modulation depth at the frequency 5oj m +2 



Hm = H (Jcj m+ 2) = A £, n 



(12) 



where A oc £l a m is a complex constant. Thus, H m is a 
measure of the complex weight £ m . For small- amplitude 
vibrations, a quantitative map of the local out-of-plane 
vibration amplitude zo can be assessed from the ratio 
of the first-order sideband Hi and the non-shifted light 
component Hq 
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where ( ) corresponds to spatial averaging over the ex- 
tent of the resonant cantilever, arrowed in fig. [2fb). In 
the same manner, a map of the local mechanical phase 
retardation ip with respect to the excitation signal can 
be assessed from Hi, and Hq 



i> - ^0 



Arg ((Hi/Hq 



(14) 



Relationships [13] and [14] are derived from eqs. [7] [J2] 
One can also make use of the third recorded band 
H-i and the property J_ 1(^0) = —Ji(<j>o) to bene- 
fit from a small SNR improvement, since we also have 

z o » it (I (^-1/^0) I) andV-V'o = Arg ((^0/^-1)). 

To perform the measurement, an optical LO consist- 
ing of the combination of three frequency components 
shifted by Acji, Aco>2, and Aujs (eqJHJ) was realized. 
Those shifts were chosen to yield non-symmetrical 
(non-opposed) algebraic modulation frequencies Slji, 
5uj2 , and 5ujs in the detector bandwidth in order to avoid 
spurious aliasing or cross-talk effects. To do so, we made 
use of phase-locked RF signals generated by frequency- 
synthesizers at the carrier frequency uoq /(2tt) = 80 MHz, 
and ujq + 5uj'i, uu + Suu^ where the shifts used to 

generate the three RF bands were Su[/ '(2tt) = —1.25 Hz, 
Suj 2 /(2tt) = 6.25 Hz, ^/(2tt) = -2.5 Hz. As sketched 
in fig. [TJ those signals were mixed and combined to drive 
a set of acousto-optic modulators designed to operate 
around 80 MHz, from which opposite diffraction orders 
were selected. The resulting modulation frequencies of 
the holograms (eqlT2]) were Suji = Suj[ + 5(jo' 2i Suj2 = 
Sujs = Suj[ — 5uj 2 . The three side band components 
H-i, Ho, and Hi, encoded at 5 Hz, -2.5 Hz, and -7.5 
Hz respectively, are visible on the magnitude of the 
discrete Fourier spectrum (eq. [T2]) reported in fig. [3] 
Spectral maps of the vibration amplitude zo (oj) were 
calculated from a set of sequential measurements at 
each frequency. We took 3000 sequences of N = 8 raw 
interferograms X at a readout rate of ujs/(2tt) = 20 Hz, 
for excitation frequencies uj/(2tt) ranging from Hz to 3 
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kHz, in 1 Hz steps (supply voltage of the piezo electric 
actuator : 10 mV). The first resonance frequencies of 
the lamellophone's cantilevers are observed between 
500 Hz and 2500 Hz (Media 1). In particular, a reso- 
nance of the first cantilever was identified around 541 Hz. 

A finer screening around this resonance was performed, 
from 248 sequences of N = 8 raw inter ferograms, for 
excitation frequencies ranging from 536 Hz to 546 Hz (the 
values of zq(uo) are reported in figHJ points). But is 
randomized from one frequency point to the next. Hence 
the phase retardation ip cannot be retrieved. To prevent 
reference phase randomization between measurements, a 
chirp of the frequency was performed. It consisted of a 
linear drift with time t of the excitation frequency from 
cji/(2tt) = 536 Hz to u f /(2tt) = 546 Hz in T = 99.2 s 

u(t) = uji + (uj f - lji) t/T (15) 

Two out of the three components of the optical LO 
(eq. E]) were chirped as well via uj(t) throughout the 
acquisition. During the time lapse T, a sequence of 
248 x 8 = 1984 raw interferograms X was steadily 
recorded. Amplitude zq{uo) and phase ip(u) maps at 
each frequency uj were calculated from the demodulation 
of N = 8 raw frame packets, reported in media 2. These 
maps were coherently averaged over the first cantilever 
according to eq. [13] and eq. [14] and reported in fig. [4] 
(blue lines). In agreement with eqHJ for non-zero vis- 
cosity, a phase shift of ~ tt is observed at the resonance 
frequency uo T ~ 541.3 Hz. The theoretical resonance lines 
in fig|4] were determined from values of ujo and £ for 
which eq|3]and eq|4]best fit the experimentally- measured 
vibration amplitude. 

In conclusion, we performed optical path length 
modulation imaging by time-averaged heterodyne holog- 
raphy in off-axis and frequency-shifting conditions. To 
do so, modulation side bands were acquired simultane- 
ously through coherent frequency-division multiplexing. 
Additionally, a linear chirp of the excitation and the 
detection frequencies allowed to map the local mechan- 
ical phase retardation of a resonant cantilever with 
respect to the excitation signal without stroboscopy. 
The proposed method enabled robust and quantitative 
imaging of out-of plane nanometric vibration amplitudes. 
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